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ABSTRACT 

v 

A 

A  new  shape  normalization  for  non-uniformly  active  isothermal  catalyst 
pellets  is  developed.  It  is  based  on  the  volume  average  of  catalyst  activity 
in  a  thin  reaction  zone  near  the  external  surface  of  the  active  catalyst 
reqion,  wherein  the  reaction  occurs  in  the  diffusion-controlled  reqime.  This 
normalization  permits  inclusion  of  those  catalyst  activity  profiles  that 
become  zero  at  the  external  surface  of  the  catalytically  active  region,  which 
were  excluded  from  previous  normalizations.  When  the  surface  activity  is 
nonzero,  results  reduce  to  those  obtained  earlier*  Ll^-2d</^~The  normalization 
provides  exactly  the  slope  of  the  effectiveness  factor  -  Hiiele  modulus 
curves,  for  all  pellet  shapes  and  activity  profiles,  in  the  limit  of 
diffusion-control.  The  development  is  initially  made  for  first-order 
reactions,  and  then  generalized  to  include  arbitrary  reaction  rate 
expressions.  Effects  of  finite  external  mass  transfer  resistances  are  also 
discussed.  4-  - . 
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SIGNIFICANCE  AND  EXPLANATION 


The  effect  of  diffusion  on  the  rate  of  chemical  reaction  in  a  catalyst 
pellet  is  evaluated  by  the  so-called  "effectiveness  factor,"  which  is  the 
ratio  of  the  actual  rate  of  reaction  and  that  in  the  ahsence  of  transport 
limitations.  For  an  isothermal  first-order  reaction,  the  effectiveness  factor 
is  a  function  of  the  pellet  geometry,  the  Thiele  modulus  (which  represents  how 
fast  the  reaction  is,  relative  to  diffusion),  and  the  Biot  number  (which 
represents  the  facility  of  external  mass  transfer  to  the  pellet,  relative  to 
diffusion).  If  the  pellet  is  uniformly  active  catalytically,  then  Aris  showed 
in  a  classic  work  [Chem.  Eng.  Sci.,  _6,  262  (1957))  that  the  length  dimension 
used  in  defining  the  Thiele  modulus  and  the  Biot  number  could  be  so  redefined, 
in  a  physical  manner,  that  the  asymptotic  behavior  of  the  effectiveness  factor 
for  large  values  of  the  Thiele  modulus  (i.e.,  the  diffusion-controlled  regime) 
was  identical  for  various  pellet  shapes.  Fortunately  with  such  a  redefini¬ 
tion,  the  effectiveness  factor  versus  Thiele  modulus  curves,  for  the  entire 
range  of  Thiele  modulus  values,  also  virtually  coincide  for  all  pellet  shapes. 
Thus,  essentially,  if  one  performs  computation  for  one  pellet  shape,  one  has 
simultaneously  obtained  the  effectiveness  factor  for  all  pellet  shapes. 

Now,  either  as  a  result  of  preparation  technigue,  which  may  well  be 
deliberate,  or  as  a  result  of  use,  catalyst  pellets  develop  non-uniform 
activity  profiles.  We  have  recently  been  attempting  (1,21  to  provide  shape- 
normalizations  in  the  above  sense,  when  the  catalyst  activity  profile  is  non- 
uniform.  Various  types  of  activity  profiles  were  previously  (1,2)  shown 
amenable  to  such  normalization.  The  present  work  brinqs  us  one  step  closer  to 
the  eventual  goal  of  including  all  possible  types  of  activity  profiles. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


ON  SHAPE  NORMALIZATION  FOR  NON-UNIFORMLY  ACTIVE  CATALYST  PELLETS  -  III 


M.  Morbidelli  and  A.  Varma 
1.  INTRODUCTION 

The  shape  normalization  for  non-uniformly  active  catalyst  pellets,  in 
which  an  isothermal  first-order  reaction  occurs,  was  examined  extensively  in 
two  previous  papers  (1,2].  Both  the  cases  of  negligible  and  finite  external 
mass  transfer  resistance  were  considered.  The  particular  case  of  catalyst 
activity  profiles  which  become  zero  in  a  continuous  manner  somewhere  within 
the  pellet  or  at  the  pellet  external  surface,  specifically  also  with 
negligible  external  mass  transfer  resistances,  was  however  previously 
excluded. 

The  purpose  of  the  present  communication  is  to  develop  a  physically  based 
shape  normalization  which  also  includes  this  particular  situation.  In  fact, 
only  the  asymptotic  slope  of  the  effectiveness  factor -Thiele  modulus  curves, 
for  larqe  values  of  the  Thiele  modulus  (i.e. ,  diffusion-controlled  regime),  is 
found  amenable  to  a  perfect  normalization.  Hie  absolute  magnitude  of  the 
effectiveness  factor,  even  for  large  values  of  the  Thiele  modulus,  is  found  to 
depend  still  on  parameters  that  characterize  the  specific  activity  profile 
under  consideration.  For  catalyst  activity  profiles  which  do  not  fit  in  the 
category  described  above,  the  normalization  developed  here  becomes  identical 
with  that  given  by  us  earlier  (2], 
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2.  ASYMPTOTIC  BEHAVIOR  AND  SHAPE  NORMALIZATION 

Consider  a  catalyst  pellet  with  an  arbitrary  given  catalyst  activity 
distribution  a(x),  in  which  an  isothermal  first-order  reaction  occurs  with 
negligible  mass  transfer  resistance  external  to  the  active  catalyst  zone.  Let 
us  now  examine  the  behavior  of  the  effectiveness  factor  for  large  values  of 
the  Thiele  modulus. 

It  is  well-known  that  in  the  diffusion-controlled  regime,  the  reaction 
occurs  near  the  surface  of  the  active  catalyst  zone,  within  a  narrow  region 
whose  width  is  proportional  to  the  reciprocal  of  the  Thiele  modulus  {31.  This 
may  also  be  illustrated  qualitatively  by  the  following  arguments. 

The  distance,  6  covered  by  a  diffusing  reactant  molecule  within  the 
pellet  can  be  evaluated  from  the  relationship 

«  ■  (Dt)'/2  (1) 

e 

where  is  the  effective  diffusion  coefficient,  and  T  is  the 

characteristic  time  available  to  the  reactant  molecule  for  diffusing,  before 
it  reacts.  A  reasonable  estimate  of  T  can  be  obtained  by  considering  a  batch 
system,  described  by  the  material  balance 

fl  -  <2» 

with  the  initial  condition:  C  «  CQ  at  t  =  0.  In  this  equation,  the 
catalyst  activity  a  is  assumed  constant  during  time.  From  the  solution  of 
eqn  (2): 

C  *  C^exp(-kat),  (3) 

the  time-constant  for  the  reaction  is 


t  -  l/(ka) 


(4) 


T 


For  the  diffusion-reaction  problem,  a  would  he  replaced  by  aj ,  the 
volume  average  of  a(x)  in  the  active  catalyst  region  where  the  reaction 
occurs i  substituting  this  in  egn  (1)  gives 

6  -  (De/ka6 )  ^ .  (5) 


Now,  introducing  the  dimensionless  width  of  the  reaction  zone  «=  6/ R, 
the  expected  relationship  with  the  Thiele  modulus  is  obtained: 


V  °e  i  V2  1 

rd  "  rI  __  J  "  J  * 


Note  that  ag  is  the  volume  average  value  of  the  catalyst  activity  in  the 
reaction  zone,  (R-d  )  <  x  <  R;  i.e. , 

JR  a(x)dV 
R-d  P 

a  .  ■  - - .  (  7 ) 

6  JR  dv 

R-d  P 

The  definition  of  the  Thiele  modulus,  eqn  (6),  thus  obtained  is  consistent 
with  that  used  previously  [1]  for  non-uniformly  active  catalyst  pellets;  the 
average  value,  a^  here  replaces  the  surface  value,  a(R)  employed  before. 

Introducing  the  dimensionless  space  coordinate  r  »  x/R,  and  notinq  from 
eqn  (6)  that  r,  ♦  0  as  $  ♦  •»,  and  therefore  dV_  *  S  dx  for  all  pellet 

0  p  X 


shapes,  eqn  (7)  reduces  to 


—  I1  a(r)dr. 
r .  . 

5  1-rfi 


The  effectiveness  factor  can  be  calculated,  by  definition,  as  the  ratio 
between  the  rate  of  reactant  consumption  in  the  actual  reaction  zone 
(1-r^  <  r  <  1),  and  the  same  quantity  in  the  case  of  fast  diffusion— in 
which  case  the  reaction  occurs  within  the'  whole  catalytically  active  region: 


(S  d )(ka,  )C 
x _ d 

v*vc. 
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where  the  subscript  a  denotes  the  asymptotic  value  of  the  effectiveness 
factor  in  the  diffusion-controlled  reqime.  a^  is,  similar  to  ,  the 
volume  average  catalyst  activity  in  the  entire  catalytically  active  region 
(which  becomes  the  reaction  zone  in  the  case  of  fast  diffusion;  i.e.,  the 
kinetic-controlled  regime).  a,)  is  therefore  given  by  eqn  (7),  by  replacing 
6  by  R,  and  this  in  dimensionless  form  reduces  to 

*1  “T^iT  i1***"*-  (10) 

Note  that  if  the  catalyst  activity  profile  a(x)  is  defined,  as  done  before 
[1,2],  by  the  ratio  between  the  local  rate  constant  k(x)  and  its  volume 
averaged  value  k,  then  by  definition,  »  1  for  all  pellet  shapes  and 
activity  profiles.  With  such  activity  profiles,  the  volume  average  reaction 
rate  constant,  in  pellets  of  different  shapes  and  different  activity  profiles, 
is  the  same.  Catalyst  activity  profiles  of  this  type  are  hereafter  called 
conservative,  since  they  are  defined  so  as  to  meet  the  condition  of  active 
catalyst  conservation. 

In  egn  (9),  the  reactant  concentration  C  in  the  reaction  zone  can  be 
estimated  as 

C  *  C  ,  (11) 

8 

since  r^  ♦  0  as  $  *  •»  and  the  reaction  zone  in  the  limit  shrinks  to  the 
active  catalyst  zone  surface.  With  this,  eqn  (9)  gives 


n 


a 


(12) 


Also,  sx/vp  *  (n+1)/R,  and  r^ 


5/R;  and  so 


If  we  now  introduce  a  normalized  ttiiele  modulus  ♦  ,  defined  as 


and  substitute  eqns  (6)  and  (14)  in  eqn  (13),  then 

n  -  1/*.  (15) 

Thus  for  different  pellet  shapes  (i.e. ,  n)  and  for  different  catalyst 
activity  profiles  a(x),  the  asymptotic  behavior  of  n  for  the  diffusion- 
controlled  regime  is  brought  together.  Note  that  for  a  conservative  activity 
profile  (i.e.,  a^  =  1),  the  expression  of  the  normalized  ttiiele  modulus 
given  by  egn  (14)  differs  from  that  used  before  [1,2]  only  in  that  a ^ ,  the 
volume  average  catalyst  activity  in  the  reaction  zone,  is  used  in  place  of  the 
catalyst  activity  at  the  outer  surface  of  the  catalytically  active  zone, 
a(R). 

From  egn  (8),  it  is  evident  that  a^is  a  function  of  the  reaction  zone 
width,  rg  and  of  the  activity  profile,  a(r).  Depending  on  the  specific 
characteristics  of  the  latter,  two  different  possibilities  can  arise,  and 
these  are  now  identified. 

Since  we  are  considering  the  case  of  $■*••,  and  therefore  r^  0,  the 
activity  profile  a(r)  may  be  approximated  by  its  Taylor  series  around 
r^  =  0  (i.e.,  r  *  1)  truncated  after  the  first-order  term: 

a(r)  -  a( 1 )  +  a'(1)rj  +  0(rg)  (16) 

where  a(1)  and  a'{1)  are  the  values  of  catalyst  activity  and  its  first 
derivative,  respectively,  evaluated  at  the  outer  surface  of  the  active 
catalyst  zone  (i.e.,  r  -  1).  Substituting  egn  (16)  in  egn  (8),  and 
performing  the  indicated  integration  qives 

a{  -  a ( 1 )  +  a'(Dr6/2  +  0(rj).  (17) 


-5- 


(18) 


If  a(1)  *  0,  then  for  r g  ♦  0  (i.e. ,  $  +  •), 

m  s  ( 1 )  • 

If,  however,  a{1)  *  0,  the  limiting  value  of  a g ,  for  r^  ♦  0,  becomes  a 
function  of  r^  and  of  the  derivative  of  the  activity  profile,  a'(l).  If 
the  latter  quantity  is  also  null,  then  a^  depends  on  the  first  nonzero 
derivative  of  a(r)  at  r  *  1. 

The  result  that  we  wish  to  stress  here  is  the  following: 

i)  if  a(1)  *  0,  then  a^  is  a  function  only  of  the  activity  profile 

itself; 

ii)  if  a(1)  *0,  then  a^  depends  also  on  r^ ,  and  therefore  through 

eqn  (6)  on  the  Thiele  modulus,  4  .  For  this  reason,  in  the  relationship 

between  n  and  $  given  by  eqns  (14)  and  (15),  the  value  of  a.  also 
a  o 

affects  the  asymptotic  exponent  of  the  Thiele  modulus,  and  not  only  the 
multiplier  as  in  the  previous  case. 


2.1.  Case  of  a(1)  *  0, 

Substituting  eqn  (18)  in  (14)  gives 


* 


Vi  f Ka(  1  )■>  V2 

S  a(  1 )  ^  D  ‘ 
x  e 


(19) 


which  is  the  same  as  that  obtained  before  (1,2).  From  eqn  (15),  n  =  1/#. 

SI 

We  do  not  pursue  this  case  further,  since  it  has  been  examined  extensively 


earlier  [1,2] 


2.2.  Case  of  a(1)  ■  0. 

In  this  ease,  5g  is  a  function  of  r ^ ,  which  depends  on  the  specific 
type  of  activity  profile  under  consideration.  We  shall  therefore  focus  on 

a(r)  »  m(1-r)®  (20) 

which  was  examined  before  (1,21.  Note  that  for  suitable  values  of  m  and 
0,  eqn  (20)  may  be  used  to  approximate  many  activity  profiles  which  are 
characterised  by  a( 1 )  ■  0. 

Substituting  eon  (20)  in  eqn  (8)  gives 

0 

a^  -  mra/{0+1)  (21) 

which  using  eqn  (6)  leads  to  an  expression  of  the  reaction  zone  width: 

r,  „„ 

0  * 


These  two  equations,  together  with  eqns  (14)  and  (15)  give 


e-2. 


n  ,  £+1  (-2L.)  V <>+2) (1SS_)- (0+D/ (0+2) 


(23) 


a1  may  be  evaluated  from  eqn  (10)  for  the  specific  a(r)  given  by  eqn  (20) 


(n+1 ) 1 
n+1 

n  (0+i) 
i=1 


(24) 


Finally,  introducing  *  R(k/D@)  ^  a  Thiele  modulus  based  on  the  average 
rate  constant  and  the  characteristic  pellet  lenoth,  eqn  (23)  may  he  rewritten 
in  the  form 


n 


a 


n+1 

_  n  (0+i) 

. -2g  (0+1 )~°  i-1 
0  g  nl 


(25) 
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where  c  *  1/10+2),  g  =  (0+1)/(0+2),  and  n  *  0,2  indicate  the  slab  and 
sphere  geometries,  respectively.  Note  that  V2<  g  <  1,  and  so  the  asymptotic 
behavior  of  the  effectiveness  factor  as  a  function  of  the  Thiele  modulus  is 
with  an  exponent  that  lies  between  -1  and  -2. 

The  exact  asymptotic  expression,  obtained  from  the  exact  analytic 

solution  (reported  in  Table  5  of  [ 2]  >  as  ♦  •  is: 

n+1 

H  (0+i) 

n  ,26) 

e  0  g  r(c)  n! 

m 

A  comparison  of  the  two  asymptotic  expressions,  given  by  eqns  (25)  and  (26), 
shows  that  the  exponents  of  the  Thiele  modulus  are  identical,  but  the 
multipliers  are  somewhat  different.  This  difference  is  examined  next,  but  it 
is  first  worth  noting  that  eqns  (25)  and  (26)  are  both  valid  for  arbitrary 
m  and  0,  whether  they  satisfy  the  conservation  condition  or  not. 


-8- 


3.  NUMERICAL  COMPARISONS 


In  order  to  determine  the  difference  between  the  two  asymptotic 
expressions  for  n ,  let  us  define  their  ratio 

n  -  n  /n  (27) 

6  d 

which,  using  eqns  (25)  and  (26)  is  given  as 

fl  «  <0+1)C(0+2)0C  ,  (28) 

a  quantity  that  is  a  function  only  of  0,  and  not  of  pellet  shape. 

It  may  now  be  seen  that  0  ♦  1  as  either  0+0  or  0  ♦  °»,  and  so  the 

two  asymptotic  expressions  match  in  both  limits.  The  value  of  £1  as  a 

function  of  the  activity  profile  parameter  0  is  shown  in  Fig.  1.  It  may  be 

observed  that  the  maximum  deviation  of  n  from  0  occurs  at  0  *  4.5, 

a  e 

where  the  relative  percentage  error  is  14.2%.  Moreover,  0  <  1  for  all 

positive  0,  which  indicates  that  h_  is  an  upper  bound  of  the  exact 

d 

asymptote  n  • 
e 

A  calculation  of  the  effectiveness  factor  n,  in  the  entire  range  of  the 

normalized  Thiele  modulus  ♦,  evaluated  by  using  the  exact  solution  (see 

Table  5  of  [2]),  is  shown  in  Fig.  2  for  the  slab  (n=0)  and  sphere  (n=2) 

geometries,  for  0=0  and  3.  As  noted  above,  it  may  be  seen  that  the 

asymptotic  slopes  of  the  various  n-$  curves  for  large  *  are  identical,  and 

egual  to  -1.  The  absolute  magnitudes  of  the  asymptotes,  however,  depend  on 

the  activity  profile  parameter  0,  but  not  on  pellet  shape.  It  may  also  be 

observed  that  the  spherical  pellet  reaches  the  asymptotic  value  at  relatively 

laraer  values  of  $.  This  is  due  to  the  presence,  in  the  exact  asymptotic 

expression  of  the  effectiveness  factor  for  the  sphere  (see  Table  5  of  [21),  of 

_2 

a  term  proportional  to  ,  in  addition  to  the  leading  term  proportional  to 


$q~~' •  Since  -2  <  -2q<  -1,  the  term  with  eventually  dominates  ♦q"? 

this  feature,  however,  delays  achievement  of  the  asymptotic  slope  -1.  Note 
also  that  this  effect  becomes  more  significant  at  larger  vlaues  of  8. 

It  is  worth  reiterating  that  when  a(1)  -  0,  the  physically  based  shape 
normalization  developed  here,  brings  only  the  asymptotic  slopes  of  the 
effectiveness  factor  -  Thiele  modulus  curves  together  for  various  pellet 
shapes  and  activity  profiles;  the  absolute  magnitude  is  not  normalized 
perfectly.  A  more  accurate  evaluation  of  the  magnitude  of  the  effectiveness 
factor  would  require  a  more  detailed  description  of  the  reactant  concentration 
profile  close  to  the  outer  surface  of  the  active  catalyst  zone.  This  means 
that  the  reaction  time,  t  would  be  evaluated  through  a  relationship  more 
accurate  than  that  given  by  eqn  (2),  which  takes  into  account  both  the 
catalyst  activity  and  the  reactant  concentration  profile  in  the  reaction 
zone.  Such  a  relationship,  however,  becomes  the  same  as  the  starting  point  of 
the  present  diffusion-reaction  problem,  and  so  the  purpose  of  an  a  priori 
shape  normalization  is  defeated. 


;  A 


4.  EFFECT  OF  EXTERNAL  MASS  TRANSFER  RESISTANCE 


Let  us  now  consider  a  pellet  in  which  the  active  catalyst  is  located  in 
an  inner  core  0  <  x  <  x^,  while  the  balance  <  x  <  R  is  an  inert 
support.  Hie  activity  profile  thus  has  the  form 


a(x)  for  0  <  x  <  x. 


for  x2  <  x  <  R 


(29) 


where  a(x)  is  an  arbitrary  qiven  distribution.  If  external  mass  transfer 
resistance  in  the  bulk  phase  is  also  taken  into  account,  the  problem  reduces 
to  the  one  considered  earlier  in  (2]. 

It  was  shown  before  (2]  that  the  effectiveness  factor  can  be  expressed  as 


n 


(30) 


where  is  the  effectiveness  factor  in  the  case  of  negligible  external  mass 

transfer  resistance  (i.e. ,  both  Bi  *  “  and  ■*“),  end  whose 

asymptotic  form  for  large  $  is  given  by  eqn  (15)  with  x2  replacing  R. 

All  the  other  quantities  are  defined  as  in  [2],  simply  replacing  catalyst 
activity  at  the  active  catalyst  zone  surface,  a(1)  by  the  average  value  in 
the  reaction  zone,  a ^ .  As  a  consequence,  besides  the  normalized  Thiele 
modulus  defined  by  eqn  (14),  the  normalized  pellet  characteristic  length  is 
also  modified,  as 

R  -  V  /(S ,a, ).  (31) 

p  2  0 

The  variable  also  enters,  through  eqn  (31),  definitions  of  the  Biot 

numbers  Bi  and  Bi.  .  Note  that  if  a(1)  *  0,  then  due  to  eqn  (18),  all 

int 

these  quantities  become  identical  to  the  correspondinq  ones  in  (2). 


The  effectiveness  factor  n ,  can  now  be  represented  as  a  function  of  its 


normalized  asymptotic  value  u,  derived  from  egns  (15)  and  (30)  as 


(ll 


« 


(32) 


for  all  pellet  shapes,  activity  profiles,  and  external  mass  transfer 

resistances,  Bi  and  Bi.  ..  All  the  h-oi  curves  reported  in  Pig,  2  of  12) 

int 

are  shown  in  Fig.  3,  including  also  the  curves  obtained  for  the  activity 

profile  (20)  with  0*0,  0.5, 1,2  and  3,  for  both  the  slab  (n*0)  and  sphere 

(n*2)  geometries,  and  the  following  set  of  values  for  Bi  and  Bi.  .  s 

int 

5,**/  **,10;  5,10.  The  consequent  enlargement  of  the  hatched  zone  is  by  the 
lower  region  in  Fig.  3.  It  may  be  observed  that  although  the  asmptotic 
behavior,  both  for  small  and  large  •  ,  is  identical  for  all  situations, 
inclusion  of  activity  profile  (20)  for  which  a(1)  *  0  causes  a  significant 
increase  in  the  differences  between  the  various  n-w  curves  for  intermediate 
values  of  u.  For  example,  the  width  of  the  hatched  region  at  w=2  increases 
from  about  33%  to  slightly  over  60%. 

Finally,  note  that  the  particular  case  when  both  Bi  and  Bi^nfc  are 
infinite,  shown  previously  in  Fig.  2,  is  excluded  in  Fig.  3. 
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5.  EXTENSION  FOR  ARBITRARY  REACTION  RATE  EXPRESSIONS 


The  procedure  leadinq  to  the  shape  normalization,  developed  here  for 
first-order  reactions,  can  also  be  extended  to  include  arbitrary  reaction  rate 
expressions,  which  exhibit  positive-order  as  C  ♦  0.  This  latter  condition  is 
satisfied  by  all  physically  reasonable  rate  expressions,  and  simply  assures 
that  the  width  of  the  reaction  zone,  near  the  surface  of  the  active  catalyst 
zone,  in  the  diffusion-controlled  regime  decreases  as  the  Thiele  modulus 
increases;  i.e. ,  r ^  ♦  0  as  4  ♦  *• 

Consider  then  a  general  reaction  rate  expression  r(c),  the  definition 
of  the  Thiele  modulus  modifies,  as  usual  [3],  to 


while  expression  (12)  for  the  effectiveness  factor  does  not  change,  since  from 
eqn  (11),  r(C)  ♦  r(Cg)  as  t  Following  the  same  steps  described  above 

for  a  first-order  reaction,  it  can  again  be  concluded  that  «  1/*  as 

In  terms  of  the  non-nor mail zed  Thiele  modulus  this  implies  that  the 

asymptotic  slope  of  the  curve  for  the  catalyst  acitvity  profile  (20)  is 

again  -2g.  This  result  is  indeed  confirmed  by  the  numerical  computations 

_  2 

reported  by  Villadsen  (41,  for  r(C)  =  XC/O+KC)  .  As  before,  the  absolute 
magnitude  of  the  asymptote  of  h  for  large  $  cannot  be  obtained  exactly 
for  this  activity  profile. 
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NOTATION 


a(x)  catalyst  activity,  k(x)/k 

a  volume  average  catalyst  acitvity 

Bi  Biot  number,  kcR/Dfi 

Bi  normalized  Biot  number,  k  R/D 

C  *  vQ  S 

Bi^nt  normalized  internal  Biot  number,  [  /  — ^  dv]~ 

c  1/(B+2)  V 2 

C  reactant  concentration 

D  effective  diffusion  coefficient 

g  (B+U/tB+2) 

k(x)  local  rate  constant 

k  volume  average  rate  eorst*.»t 

kc  external  mass  transfer  coefficient 

m  parameter  characterizing  activity  profile 

n  integer  characteristic  of  pellet  geometry*  0  for  infinite  slab, 

1  for  infinite  cylinder,  2  for  sphere 

r ^  dimensionless  width  of  the  reaction  zone 

R  characteristic  pellet  dimension,  half-thiehness  (n*0),  radius 

(n-1,2) 

R  normalized  pellet  dimension, 

S  pellet  surface  area  at  location  x 

S2  pellet  surface  area  at  x  -  x2 

Sx  pellet  external  surface  area 

Vp  pellet  volume 

x  distance  from  center  of  pellet*  regions  0  <  x  <  x1  and  x  *  x  <  R 

are  inert,  while  region  x1  <  x  <  x2  is  catalytically  active. 

y  normalized  distance  from  center,  x/R 

y0  R/R 

-18- 


_ _ ‘  ~  »  ,  .. 

*s— "  I*  ■  >  nr  iBm  r  i  i  "  l—jj 


-Aik 


'3  -  ■  - 


Greek  Symbols 

6  parameter  character! zinq  the  activity  profile 

5  width  of  the  reaction  zone 

T (  )  Gamma  function 
n  effectiveness  factor 

hD  effectiveness  factor  for  the  Dirichlet  problem 

V, 

♦  Thiele  modulus,  R{ka^/De)  'z 

A  Thiele  modulus  based  on  volume  averaqe  rate  constant,  R(k/D  )  ^2 

0  © 

•  normalized  Thiele  modulus,  defined  by  eqn  (14) 

t  characteristic  time  for  diffusion 

w  parameter  defined  by  eqn  (32) 


n 


Vn 


a 


Subscripts 

a  asymptotic  value 

e  asymptotic  value  of  exact  analytic  expression 

s  surface  of  catalytically  active  zone 


6 


reaction  zone 


"T 
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